Magnetic impurities on the surface of a topological insulator 
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The surface states of a topological insulator are described by an emergent relativistic massless 
Dirac equation in 2 + 1 dimensions. In contrast to graphene, there is an odd number of Dirac points, 
and the electron spin is directly coupled to the momentum. We show that a magnetic impurity 
opens up a local gap and suppresses the local density of states. Furthermore, the Dirac electronic 
states mediate an RKKY interaction among the magnetic impurities which is always ferromagnetic, 
whenever the chemical potential lies near the Dirac point. These effects can be directly measured 
in STM experiments. We also study the case of quenched disorder through a renormalization group 
analysis. 
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Following the recent theoretical prediction and the ex- 
perimental observation of the quantum spin Hall insu- 
lator state in two dimensions [l|, 0, H, Hl^ the concept of 
a topological insulator (TI) in three dimensions has at- 
tracted a lot of interest^, @, 0, S]- The electronic exci- 
tation spectrum of a time reversal invariant TI is fully 
gapped in the bulk, but there are gapless surface states 
described by the 2 + 1 dimensional relativistic Dirac equa- 
tion with an odd number of Dirac points. This property 
makes the surface system an extremely unusual (2 + 1)- 
d system, just like the 1 + 1 dimensional "helical" edge 
states of quantum spin Hall insulators [^,[l3]- In fact, one 
can prove a general no-go theorem, which states that a 
two dimensional time reversal invariant lattice model can 
not have an odd number of Dirac points [9]. For exam- 
ple, the familiar graphene model on a honeycomb lattice 
has four Dirac points [ll|, |l2[. The surface states of a TI 
evade this no-go theorem since they describe the bound- 
ary of a three-dimensional lattice model, and a pair of 
Dirac points can be separated onto the two opposite sur- 
faces. The Dirac points of a TI are thus stable and robust. 
They can not be destroyed by any time reversal invari- 
ant perturbations. In contrast, since the Dirac points 
of a two-dimensional lattice model occur in pairs, they 
can be pairwise annihilated by small perturbations. For 
example, a sublattice distortion in graphene can remove 
the Dirac points entirely. 

Therefore, the surface states of a TI offer a unique 
platform to investigate the physics of robust Dirac points. 
In Ref [S,, J3.], it was pointed out that a time reversal 
breaking perturbation on the surface is the most natural 
way to reveal the topological properties of Dirac points. 
For this reason, we investigate the effects of magnetic 
impurities on the surface states of a TI. We consider the 
simplest case of a single Dirac point, described by the 
low-energy effective Hamiltonian 

(1) 



where the z-direction is perpendicular to the surface and 
Vf is the Fermi velocity. At first sight, this is exactly the 
2D Dirac Hamiltonian at one nodal point of graphene 
which has been used successfully to describe its low- 
energv physics 0,111. However, there is one important 
difference between these two cases. For graphene, the 
two components of the Dirac Hamiltonian describe the 
two sublattices or pseudo-spin degrees of freedom, while 
in the case of a TI, the two components describe the real 
electron spin, and are related to each other by time rever- 
sal. Therefore, we expect the coupling between magnetic 
impurities and electron spin to take the form 

He. = m,+Hl = ^ J.s.(r)5,(r)+J||(s,5,+s,5,)(r), 

r 

(2) 

where 5'i(r) is the spin of a magnetic impurity located at 
r, Si{r) = ip\r)a'^ip{r) is the spin of the surface electrons 
and J^, J|| are the coupling parameters. The Hamiltoni- 
ans ([1]) and ([2]) together describe the problem of magnetic 
impurities on the surface of a TI, which is the starting 
point of this paper. This phenomenological Hamiltonian 
can be derived rigorously from more realistic models. For 
example, strained HgTe, which is expected to be a TI, can 
be described by a realistic Kane model [1, [ij • A straight- 
forward derivation, following the similar steps outlined 
in Ref. [8], gives Hvf ^ 3.4 eV ■ A, ^ -20 meV and 
J|| ^ -40 meV. 

In graphene, the chemical potential is automatically 
tuned to the Dirac point for a half-filled system. In a 
TI the chemical potential in general needs to be tuned 
to the Dirac point by a metallic gate, which screens the 
Coulomb interaction between the surface electrons to an 
irrelevant local four fermion interaction. In our current 
work the Coulomb interaction will be ignored through- 
out, and our results are applicable to length scales larger 
than the screening length. 

Single magnetic impurity.- Let us start by studying 
the effect of a single magnetic impurity on the surface 
states. For simplicity, we assume that the impurity is 
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located at the origin and treat it as a classical spin. 
Under such a mean-field approximation, the exchange 
Hamiltonian is written as Hex — Xlr ^«(^)"^«(^) ^^^^ 
Mi(y) = Ji {Si)S{r). It should be noticed that the S- 
function form of the exchange interaction is only correct 
in the long- wavelength limit. For the description of local 
physics the J-function potential should be regularized as 
a finite range potential. In the following, we take a square 
well regularization for the exchange interaction and study 
the case with the impurity spin polarized in the z direc- 
tion. In this case the exchange Hamiltonian is reduced 
to 

r 

where M(r) = Mo6(ro — r) with 6(ro — r) the step 
function and Mq = Jz {Sz)- determines the range of 
the exchange interaction. This problem has azimuthal 
symmetry and can be solved analytically with Bessel 
functions [15]. From the analytic solution we find that the 
wavefunction for r < ro decays for energies |^| < |Mo|, 
and oscillates for energies \E\ > \Mo\. To study ob- 
servable consequences of this impurity effect, we cal- 
culate the local density of states (LDOS) defined by 
Po{r,E) = -^^{Tr[G^(r,r, £;)]}, where G^(r,r^£;) 
is the retarded single-particle Green's function. As shown 
in Fig. [T](a) and (b), the LDOS is suppressed in the en- 
ergy range \E\ < |Mo| and spatial range r < tq. For 
r > To the LDOS converges quickly to the impurity-free 

\E\ 

value ^ '2 2 • Such a LDOS gap induced by a magnetic 

impurity can be observed by STM experiments, and de- 
fine a sharp criterion to distinguish the TI surface from 
other two-dimensional systems with an even number of 
Dirac cones, such as graphene. In graphene, the two com- 
ponents of the Dirac Hamiltonian represent the pseudo- 
spin degree of freedom of the two sub-lattices, which 
don't couple to magnetic impurities directly. Therefore, 
no suppression of LDOS win be observed 16 . 

Another interesting physical quantity is the spin 
LDOS, defined by pi (r, E) = -^<^{Tr[G^ (r, r, E) a^]}, 
i = X, 7/, z. Experimentally, the spin LDOS can be 
measured by the recently developed spin-resolved STM 
technique In order to calculate the spin LDOS in- 
duced by a magnetic impurity with arbitrary polariza- 
tion, we can not use the analytic solution discussed above 
since rotation symmetry is broken. Therefore, we use 
the T- matrix formalism [18j to calculate the Green's func- 
tion (r, r\ E). The calculated distribution of the spin 
LDOS Pi (r, E) at a given energy is shown in Fig. [T] (c) 
for the out-of-plane and (d) in-plane magnetization of the 
impurity. As seen from Fig. [I](c), the z-direction mag- 
netization induces not only a z-direction spin LDOS, but 
also an in-plane spin LDOS. This is a direct consequence 
of the spin-orbit coupling of the surface states. In fact, 
the Dirac Hamiltonian ([1]) can be regarded as an elec- 
tron spin coupled to a momentum-dependent effective 
magnetic field Bg// = hvf'k. As an electron propagates, 
its spin processes around Bg//. Since Bg// is parallel to 
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FIG. 1: (a) Charge local density of states (LDOS) as a func- 
tion of the distance r from the magnetic impurity for electron 
energies ^ = 30 meV (black line) and E = 70 meV (blue 
line), (b) Charge LDOS as a function of electron energy E 
at positions r = (blue line) and r = 20 nm (red line). Here 
we assume a magnetic impurity strength Mo = 50 meV and 
a coupling range ro = 13 nm. The spin LDOS is plotted as 
a function of position ai E = 10 meV for magnetization of 
the magnetic impurity placed at the origin (0,0) in (c) the 
z-direction and (d) the y-direction. Here the arrow indicates 
the in-plane spin LDOS and the color shows the 2;-direction 
spin LDOS. 



the propagation direction of the electron, we immediately 
realize that the spin of the electron will always precess 
in a plane perpendicular to the direction of propagation. 
For example, when an electron with spin polarized in the 
z-direction is moving towards the x direction, the spin is 
precessing in the yz plane. Consequently, the spin LDOS 
vector ps{r,E) = {px, py^ pz){r, E) is canted towards the 
^/-direction when r is along the positive x-axis. A similar 
analysis can apply to other directions, from which we can 
understand the spin LDOS distribution in Fig. [T](c). For 
an in-plane impurity spin polarized in the ^-direction, the 
precession picture remains valid. For example, no preces- 
sion occurs along the ?/-axis since the spin is precessing 
in the xz plane, in agreement with Fig. [1] (d). The spin 
LDOS pattern in Fig. [1] (c) and (d) can also be under- 
stood by symmetry. When the impurity spin is polarized 
in the z-direction, the system is invariant under the ro- 
tations about z-axis, and so is the spin LDOS. For an 
in-plane impurity spin, rotation symmetry is broken, but 
there is a discrete symmetry defined by a 7r-rotation along 
the z-axis combined with a time-reversal transformation. 
Such a residual symmetry is preserved in the pattern of 
Fig. [11(d). 

Besides the patterns discussed above which are due to 
spin precession, another important feature of the spin 
LDOS is the longitudinal decay ex 1/r^. As a result 



of integrating over energies below the Fermi energy, the 
local spin polarization (s(r)) will behave as Such 
a 1 /r^ power law is a direct consequence of the fact that 
the spatial scaling dimension of the spin density in the 
free Dirac fermion theory is 2. 

Random magnetic impurities.- In the following we 
will focus on the behavior of the system with randomly 
distributed magnetic impurities on the surface S(r) = 
J]]. Si(5(r — R^), where are the positions of magnetic 
impurities. As every magnetic impurity will open a local 
gap in its vicinity, we may expect the system to be gapped 
everywhere, at least at the mean field level. However, this 
is not necessarily true if the magnetization of the mag- 
netic impurities is nonuniform. To see this, we consider 
again the mean-field form of the exchange Hamiltonian 
(j3|) with a magnetization domain wall along the ?/-axis 
at X = 0, which is given by M{x) > (M(x) < 0) 
for X > (x < 0). Solving the Schrodinger equation 
directly on the domain wall, we obtain gapless chiral 
fermion modes along the domain wall with wave func- 
tion t/j ~ {l^i)^ ex.-p[ikyy — ^^^^ dx] and energy dis- 
persion E = hvfky. Thus, this system is in fact not 
totally gapped but has gapless modes. Compared with 
the fully gapped system, the appearance of such gapless 
modes will cost more energy. Therefore, heuristically we 
expect the system not to favor any magnetic domain wall, 
which indicates that magnetic impurities should be fer- 
romagnetically coupled. 

Keeping such a heuristic picture in mind, we now study 
the coupling between two magnetic impurities microscop- 
ically. The itinerant electrons can mediate a spin inter- 
action between two magnetic impurities, known as the 
Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction. 
Such a coupling can be obtained by integrating out the 
fermions in the Hamiltonians ([1]) and ([2]), which results in 
the form H^n = Ei,,=.,,,. - r'|)5u(r)52,(r') for 

any two magnetic impurities Si and 5*2. The coupling 
constant ^ij{R) is a function of R = |r — r^| and can be 
extracted from standard second-order perturbation the- 
ory. For example, the z-direction coupling constant has 

the form ^,,{R) = ^M.{F^{kFR) ^ F^ikpR)), where 
/cf is the Fermi momentum, ao is the lattice constant and 

— {-\-)Ji{x)Ji{x^)]. Jn{x) is the Bessel function and 
Xc = kcR with kc a large momentum cutoff. The os- 
cillating part of the RKKY interaction is determined 
by + F_ and the decaying part is proportional to 
1/R^. Such dependence is related to the depen- 
dence of the local spin polarization induced by a sin- 
gle magnetic impurity discussed earlier. Being a conse- 
quence of the Dirac Hamiltonian, similar behavior has 
also been found in graphene[19, 20]. Moreover, and 
F- give the intra-band and inter-band contributions re- 
spectively. The novel property of Dirac fermions appears 
when the chemical potential is close to the Dirac point. 
The oscillation period of the RKKY interaction being 
determined by the Fermi wavelength Xp = l/kp^ the os- 




FIG. 2: RKKY interaction versus the distance R between 
two magnetic impurities. Fermi momentum /cf is chosen to 
be 0.5/ao for the blue line, 1.0/ao for the green line and 1.5/ao 
for the red line, where ao is the lattice constant. 



dilation becomes weaker as /c^ ^ 0, as shown in Fig. 
[21 Eventually the two magnetic impurities become fer- 
romagnetically coupled when Af is much larger than the 
average distance between them. The RKKY interaction 
mediated by the surface states of a TI is quite differ- 
ent from that between magnetic impurities in graphene. 
The tendency towards ferromagnetic [2l|,[22[ or antiferro- 
magnetic correlations [l^ HO, [23] in graphene is a topic 
of current debate. Conside ring only the noninter acting 
lattice model for graphene [Ull^, it is found that the 
RKKY interaction is ferromagnetic for local moments 
within equivalent sublattices but antiferromagnetic for 
opposite sublattices. Such different behavior is again a 
consequence of the different physical meaning of the two 
components of the Dirac spinor in the two systems. 

Due to the ferromagnetic RKKY interaction discussed 
above, we expect ferromagnetic order to appear in this 
system when the chemical potential is near the Dirac 
point. At the mean field level, from the H^^ term 
in ([2]) we know that Jz{sz) acts as an effective mag- 
netic field to magnetize the magnetic impurities. At the 
same time, Jz{Sz) acts as the effective magnetic field 
to polarize the electron spin. The behavior of the crit- 
ical temperature Tc for the ferromagnetism can be ex- 
tracted from a standard calculation [24], and is given by 

ksTc = ^^^!^^j^Jz{Ec - Ef), where E, is a cut- 
off energy, Ef is the Fermi energy, and Sq is the sat- 
uration spin value of each magnetic impurity. Setting 
y 0.1, = Q A, Ec = 0.1 eV and 5/2 for 

Mn^+, Tc is estimated to be of the order of ~ 0.1 
K. Interestingly, such a ferromagnetically ordered sur- 
face state also carries a half quantized Hall conductance 
gh = ±e^/2/i, which can be understood as the parity 
anomaly of massless Dirac fermions [25*, ^26^, and is a di- 
rect manifestation of the nontrivial topology of the bulk 
system. To go beyond the mean-field picture studied 
above, we now switch gears to a renormalization group 
(RG) approach with more generic quenched disorder. 
We consider three types of random potentials: random 
charge potential coupled to V^^V^, random XY- moment 
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coupled to ip^a^^y^ip, and random Z-moment coupled to 
tp'^a^tp. The La grangian density including all the random 
potentials is £ = '''(Sr+^'^/cr-^ (r )?/;''' (r)?/;(r) — 

Vfg2B2{r)jtlj^r)a^tlj{r) vjgsBs{r)tlj^r)a^tlj{r) , where 
j = x^y. Ba{v) represents three types of random poten- 
tials with correlation {B{y)B{y')) = d'^{T - r')/(27r). gl 
is proportional to the standard deviation of each type of 
random potential distribution, which from naive power- 
counting is a marginal perturbation. Since the random 
potential has long range temporal correlation but short 
range spatial correlation, Lorentz invariance of the bare 
Dirac Lagrangian is lost, and the Fermi velocity Vf flows 
under RG. In the calculation we assume that the stan- 
dard deviation ~ g^ is small. After integrating out 
degrees of freedom between momentum cutoff A and A, 
the leading order RG equations read 
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is the wavefunction renormalization and Z^o, Zyxy 
and Zyz are vertex corrections to V^^V^, V^V^'^^^V^ and 
ip^a^ip respectively. The complete solution of the RG 
equations (|4]) is complicated, but the solutions in different 
limits can be obtained straightforwardly. For instance, if 
9i = 93 = 0, the solution of the RG equations reads: 

^ = 1 + 1^, ^H = |^, A[^V]=2, 
A[iP^a'^P] = 2, A[^/;V^(^V] = 2 + (5) 

The coupling constant g2 does not flow under RG be- 
cause all the corrections can be absorbed into Vf. The 
dynamical scaling dimension z is changed because the 
Fermi velocity acquires a nonzero scaling dimension un- 
der RG, leading to the following scaling of the LDOS: 
p{uj) ~ [13]. Using the scaling dimensions of the 



fermion bilinears listed in Eq. [5l the long distance spin 
distribution around an isolated magnetic impurity with 
in-plane moment is given by /^^^^^(^^(r) : 5^(2^) (r) = 
V^tcr^(2/)^(r) ^ s+g'i/.^ ' Thus the spin LDOS pattern 
remains similar to that shown in Fig. [1] (c) and (d), but 
decays with a different power law. 

Quenched disorder has been studied ingra^hene, with 
N = 4 flavors of Dirac fermions [H, [23, llo]- In this 
case, ripples of a graphene sheet, interpreted as a random 
"gauge potential" B2{x)jilj^a^ilj^ have attracted most of 
the attention. In contrast with our results, the 1/N ex- 
pansion is usually taken when the RG equations are de- 
rived for graphene. The competition between random 
potentials and the Coulomb interaction have also been 
studied in graphene. In this case, the RG equations lead 
to various nontrivial fixed points |3l| . 

In conclusion, we have investigated the effects of mag- 
netic impurities on the surface states of a TI. A magnetic 
impurity breaks time reversal symmetry and suppresses 
the low energy LDOS locally. The surface states medi- 
ate a coupling between the magnetic impurities which 
is always ferromagnetic when the chemical potential lies 
close to the Dirac point. Therefore, we expect that a fi- 
nite concentration of magnetic impurities would give rise 
to a ferromagnetic ground state on the surface. This 
mechanism provides a physical realization of the novel 
topological magneto-electric effect discussed in Ref. [H, 
which requires breaking of time reversal symmetry on 
the surface of a TI. We also investigated the effect of 
quenched impurities on the surface states and presented 
the RG equations governing the flow of the coupling con- 
stants. The distinct signatures of magnetic impurities 
on the surface states of a TI discussed in this work can 
be readily observed in STM experiments, possibly on the 
surface of BixSbi_x [32| or strained HgTe. 
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